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We study the heavy-heavy-light quark (QQq) system in a non-relativistic potential model, and 
investigate the quark motional effect on the inter-two-quark potential in baryons. We adopt the 
Hamiltonian with the static three-quark potential which is obtained by the first-principle calculation 
of lattice QCD, rather than the two-body force in ordinary quark models. Using the renormalization- 
group inspired variational method in discretized space, we calculate the ground-state energy of QQq 
systems and the light-quark spatial distribution. We find that the effective string tension between 
the two heavy quarks is reduced compared to the static three-quark case. This reduction of the 
effective string tension originates from the geometrical difference between the inter-quark distance 
and the flux-tube length, and is conjectured to be a general property for baryons. 
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I. INTRODUCTION 

In hadron physics, the intcr-quark interaction is one of 
essential properties reflecting the non-perturbative gluon 
dynamics based on the SU(3) gauge symmetry. Par- 
ticularly the quark confinement is a nontrivial, non- 
perturbative, and unsolved problem even in static quark 
cases. The interaction between finite-mass quarks would 
include more complicated effects. Additionally, in con- 
trast to mesons, quarks in baryons have various kinds of 
the motion, configuration, and so on, and such a quark 
degree of freedom would be important for the inter-quark 
interaction. The inter-two-quark potential in baryons can 
be influenced by such nontrivial effects of the other quark. 

To investigate such finite-mass-quark "motional" ef- 
fects on the effective inter-two-quark potential, we study 
the heavy-heavy-light quark (QQq) system in this work. 
We treat the two heavy quarks in the QQq system as 
infinitely heavy, keeping the light-quark mass finite. We 
artificially change the inter-heavy-quark distance R, and 
calculate the energy of the QQq system as a function of 
R, which we call the QQq potential. This QQq potential 
includes not only the gluonic effect but also the nontriv- 
ial light-quark effect, and would differ from the static 
three-quark (3Q) potential. 

From our previous lattice QCD study the QQq 
potential is almost described as a linear and one-gluon- 
exchange Coulomb potential, like the quark-antiquark 
(QQ) potential. However, the effective string tension be- 
tween two heavy quarks is reduced compared to the static 
QQ or 3Q case, 0.89 GeV/fm. This reduction would 
be the result from nontrivial light-quark effects. Lattice 
QCD is a powerful tool for calculating hadron masses, 
potentials, and several QCD properties, but has a diffi- 
culty to clarify the light-quark wave function. Instead, 
in this paper, we employ a simple potential model, and 
calculate the QQq potential and the spatial distribution 
of the light-quark wave function, in order to explain the 
reduction mechanism of the effective string tension. 

In most of ordinary quark models and potential mod- 
els, the quark confining force in baryons is the simple 



two-body force 0, H Q. However, it is recently found 
that the confinement potential in baryons is expressed 
by the Y-type flux tube picture i, i, M H S Efl] • For 
more accurate calculations, quark models and potential 
models need to treat the confinement force in baryons 
as the three-body force. In particular, since we are now 
interested in the light-quark effects on the inter-quark 
potential, the accurate treatment for the three-quark in- 
teraction is important. Hence, we adopt the three-quark 
potential from the lattice QCD result, which is the non- 
perturbative first-principle calculation. 

As an example of realistic QQq baryons, the doubly 
charmed baryon is experimentally observed. In 2002, 
the SELEX Collaboration at Fermilab observed (dec) 
through a weak decay S+ — > A+if~7r+ [Til ] . They also 
confirmed another decay process S+ — ► pD + K~ [l2l ]. 
In their experiments, the mass is measured about 3519 
MeV. Doubly charmed baryons are also theoretically in- 
vestigated in lattice QCD [HI], potential model 0, flij . 
and other approaches [1 51 ] . 

We already know analogous systems to this QQq 
baryon, for example the H^ ion in molecular physics, 
which is also the bound state of two heavy and one 
light particles. However, the interaction in the ground- 
state QQq system is attractive, while the Hj ion includes 
Coulomb repulsive force between the two protons. Addi- 
tionally, as stated above, the confining force in the QQq 
baryon is a purely three-body force, not a sum of two- 
body forces, and has the characteristic geometrical struc- 
ture. Then the QQq potential is expected to include non- 
trivial effects from this feature. 

In this paper, we study finite-mass-quark effects on the 
inter-two-quark potential in baryons through the inves- 
tigation of QQq systems with a non-relativistic potential 
model. The article is organized as follows. In Sec. II, 
we introduce the Hamiltonian with the confinement po- 
tential obtained from recent lattice QCD calculations, 
and present the formalism of the renormalization-group 
inspired variational calculation in a discretized space to 
solve the Schrodinger equation. In Sec. Ill, we show the 
numerical results of the QQq potential, and discuss the 
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finite-mass-quark effect for the reduction of the effective 
string tension between the two heavy quarks in QQq sys- 
tems. Sec. IV is devoted to the summary and the con- 
clusion. 




II. FORMALISM 
A. Hamiltonian 

The Hamiltonian of three quarks in baryons is con- 
structed from kinetic terms and a three-quark interaction 
term as 

3 

# = ^T i + V(f 1 ,f 2 ,f 3 ). (1) 

In the QQq system, we treat two heavy quarks as in- 
finitely heavy particles, and one light quark as a non- 
relativistic constituent quark with the constituent mass 
M q . The non-relativistic quark model is one of the suc- 
cessful models to describe baryons even for the light- 
quark sector, and have been used for the study of 
baryons by many theoretical physicists even at present 
[l6j . Apart from irrelevant constants, the Hamiltonian is 
simplified as 



H = M„ 



_______ 

2M q dr% 



+ V(fi,r 2 ,f 3 ), 



(2) 



where the subscripts 1, 2, and 3 mean the two heavy 
quarks and the light quark, respectively. Although we 
adopt the non-relativistic formalism, we still call this 
finite-mass quark the "light" quark in this paper. 

As the three-quark interaction V{f\, r_, r 3 ), to treat 
the light-quark effect as precisely as possible, we adopt 
the lattice QCD result of the static 3Q potential [f|. This 
3Q potential includes the confining potential as the three- 
body force, instead of the simple sum of the two-body 
force in ordinary quark models. The static 3Q potential 
obtained by quenched lattice QCD is 



V(fi,f 2 , f 3 ) = a 3Q L min - ^2 — ~ 



O. 



3Qi 



°3Q 



where r,-, = 



i<j 

0.89 GeV/fm, 



.4 



3Q 



0.13, 



(3) 
(4) 

and these values in Eq.((4]) are re- 
lated to the QQ case as <t 3 q ~ °qq an( i ^3Q — 
P|. Since we are not interested in a constant shift of the 
energy in the non-relativistic formalism, we set C 3 q = 0. 
The symbol L m { n is the length of the color flux tube min- 
imally connecting the three quarks, which is described as 
follows. When all the angles of the 3Q triangle is less 
than 27r/3, the Y-type flux tube is formed and 

1 

+ \/3( r i2 + r i3 + r 23 )(-r 12 + r 13 + r 23 ) 

n 1/2 

x V O12 - r 13 + r 23 )(r 12 + r 13 - r 23 ) ■ (5) 
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FIG. 1: The schematic figure of the renormalization group 
(RG) inspired variational calculation. The finer-mesh calcu- 
lation is done with the initial condition constructed from the 
rougher-mesh result. 



When one of the angles of the 3Q triangle exceeds 27r/3, 



imin = r 12 + r 13 + r 23 - max(ri 2 , r 13l r 23 ). 



(6) 



Once the heavy-quark coordinates ri and r 2 are fixed, the 
interaction depends only on the light-quark coordinate 
f 3 . We can calculate the ground-state light-quark wave 
function ^(^3) with the variational principle of the energy 
of the system 



E(R) 



J d 3 r 3 r{r 3 W{r 3 ) 
f <Pr 3 \i>(r 3 )\ 2 



(7) 



We determine the ground-state QQq potential VQQ q (R) 
by minimizing E(R). 

Here, we comment on the three-quark interaction 
V(fi,f 2 ,r 3 ) with finite-mass quarks. The finite-mass 
effect changes the simple Coulomb interaction to the 
Fermi-Breit interaction [14| , which includes spin-spin and 
spin-orbit interactions. These relativistic corrections are 
suppressed by the inverse of the quark mass. In the QQq 
system with infinitely heavy quarks, most of them give 
zero contributions as 1/Mq — > 0, and the interaction re- 
mains the simple form. However, the finite-mass effect 
on the quark confinement potential is unknown. Here, 
for simplicity, we assume that V(ri,r 2 ,r 3 ) in the QQq 
system can be written with the static one. Equivalently, 
the finite-mass effect is assumed to be taken only via the 
light-quark wave function spreading. 



The renormalization-group (RG) inspired 
variational calculation 



We exactly solve the energy variational problem in dis- 
cretized space. We take a cylindrical coordinate (/?, 9, z), 
and locate the two heavy quarks on (0, 0, R/2) and 
(0,0, — R/2). The ground-state light-quark wave func- 
tion is mirror symmetric to the z = plane and rota- 
tional symmetric around the z-axis. Thus, we have only 
to calculate on the two-dimensional plane (p, z) (p > 0, 
z > 0). We discretize the space with an "isotropic" mesh 
as Ap = Az, and vary the light-quark wave function on 
each sites to minimize E(R). This is equivalent to solving 
the Schrodinger equation exactly, and we have no Ansatz 
about the functional form of the light-quark wave func- 
tion. 
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FIG. 2: The light-quark spatial distribution \ip(ra)\ 2 for M q = 330 MeV with R = 0.06 fm (left), R = 1.2 fm (center), and 
R — 2.4 fm (right). The brighter region has higher probability, and the black circles denote the positions of the heavy quarks. 
The figure is a part of the whole volume, and the actual calculation is performed in enough large volume. 



TABLE I: Numerical results for the different light-quark masses M q and the different mesh sizes Az(= Ap). The calculation 
listed here is done with the 128 x 128-mesh. The omitted length unit is fm. 



M q [GeV] 


Az 


y^)(iJ=0.06) 


^)(R=1.2) 


y^y(i?=0.06) 


^)(R=1.2) 


(Lmin) (#=0.06) 


(Lmin) (R— 


0.33 


0.08 


0.41 


0.42 


0.43 


0.53 


0.67 


1.61 




0.10 


0.42 


0.41 


0.45 


0.53 


0.67 


1.60 




0.12 


0.40 


0.40 


0.46 


0.55 


0.68 


1.60 


0.50 


0.08 


0.36 


0.36 


0.39 


0.47 


0.58 


1.53 




0.10 


0.35 


0.35 


0.39 


0.49 


0.59 


1.53 


1.0 


0.05 


0.27 


0.29 


0.30 


0.42 


0.46 


1.44 




0.10 


0.27 


0.28 


0.32 


0.42 


0.47 


1.43 


2.0 


0.05 


0.20 


0.23 


0.23 


0.37 


0.35 


1.36 




0.08 


0.20 


0.22 


0.24 


0.37 


0.36 


1.36 



We adopt the variational calculation inspired by the 
renormalization group (RG) method. The schematic pro- 
cedure is shown in Fig. [1] and its concrete process is as 
follows. First, we start with a 2x2 mesh and the spacings 
ApW and Az {1 \ We minimize E(R), and then obtain 
the 2x2-solution ^>W(Z,m). Here the site (l,m) corre- 
sponds to (p,z) = (2ApW,mAzW) (l,m g Z). Next, 
we turn to a 4x4 mesh with the twice finer mesh size, 
starting from the initial condition ip^ '(I ,m) from the 
2x2-mesh solution. The 2 n+1 x 2" +1 -mesh initial condi- 
tion is set from the 2™ x 2™-mesh solution, as 



4™ +1) (2Z,2m- 1) 

= ~{^(l,m) + ^ n \l,m-l)}, (12) 
^™ +1) (2Z-l,2m- 1) 

= ~{<0 (n) (Z,m) + V (n) (/- l,m- 1)}. (13) 



A P {n+1) = \^P (n) , (8) 



( n+ -n 1 /„•) We repeat this procedure N times, finally obtain the 

Az (n = > ( 9 ) 2 N x 2 JV -mesh solution with the spacings Ap = ApW 

(21 2m) — il)( n ^(l m) (10) &n( ^ ^ Z ^ ^-^^ ■ With this RG inspired variational cal- 

k ' /' v ) culation, the solution is expected to converge rapidly to 

ipQ +1 ' (21 — 1, 2m) the absolute minimum. To estimate the discretization er- 

l . . ror and the finite-volume effect, we calculate with several 

= ttIV' (l,m)+ip (I — l,m)}, (11) mesh sizes and mesh numbers. 
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FIG. 3: The QQq potential VQQ q (R) from the variational 
calculation. The different four symbols denote the result with 
the different four light-quark masses M q = 0.33, 0.50, f .0, 2.0 
GeV. The solid curves are the best-fit functions of Eq. (|f4p . 



III. RESULTS 
A. The QQq potential 

Three examples of the obtained light-quark spatial dis- 
tribution are shown in Fig.O When R becomes large, the 
light-quark distribution is broadened in the z direction. 
To estim ate t he size of th e ligh t-quark spreading, we cal- 
culate y/ (z 2 ) and yj (x 2 ) = y (p 2 )/2 from the obtained 
wave function. Some typical examples of these values are 
shown in Table |TJ \J (x 2 ) is almost i?-independent, and 
is about 0.4 fm for M 9 =330 MeV. 

We fit the QQq potential with 

V QQq (R)=a eS R-^+C eS (14) 

as an analogy of the QQ potential. The "eff" means effec- 
tive values including the light-quark effect. The best-fit 
parameters are summarized in Table [TT1 and the obtained 
QQq potential is shown in Fig. [3J This fitting function 
is significantly suitable for Vqq 9 . The effective Coulomb 
coefficient A e g is almost the same value as A^q. The ef- 
fective string tension a e g is reduced about 10-20% com- 
pared to the string tension <jzq in the 3Q potential, 

a cS < a 3Q ~ <7qq ~ 0.89 GeV/fm. (15) 

This result means that the inter-two-quark confining 
force in baryons is reduced due to the light-quark ex- 
istence. 

To estimate the discretization error and the finite- 
volume effect, the results with different mesh numbers 



TABLE II: The best-fit parameters <r e a, A c s, and C c s in 
Eq. (TJU with the fit range of R < 1.2 fm. The 128x128- 
mesh is used. This C c s differs from C e a of the lattice result 
in Ref. Q. 



M q [GeV] 


Az [fm] 


o-ea [GeV/fm] 


A eS 


CW [GeV] 


0.33 


0.08 


0.74 


0.12 


1.06 




0.10 


0.73 


0.12 


1.07 




0.12 


0.73 


0.12 


1.07 


0.50 


0.08 


0.75 


0.12 


1.11 




0.10 


0.75 


0.12 


1.10 


1.0 


0.05 


0.80 


0.12 


1.41 




0.10 


0.78 


0.12 


1.42 


2.0 


0.05 


0.84 


0.12 


2.23 




0.08 


0.83 


0.11 


2.24 


TABLE III: The different-mesh results with M, 
o"eff, A e s, and C e ff are the best-fit parameters in 
the fit range of R < 1.2 fm. 


, = 330 MeV. 
Eq. CO} with 


Mesh 


Az [fm] 


(Jcft [GeV/fm] 


A eS 


C cff [GeV] 


64x64 


0.12 


0.72 


0.12 


1.07 




0.15 


0.72 


0.12 


1.08 


128x128 


0.08 


0.74 


0.12 


1.06 




0.10 


0.73 


0.12 


1.07 




0.12 


0.73 


0.12 


1.07 


256x256 


0.05 


0.74 


0.12 


1.08 



and mesh sizes are shown in Table IIII1 These results 
are different only within a few percent. We can see that 
our variational calculation is performed in enough large 
volume and enough fine mesh. 

In our previous lattice QCD study [ij, we investigated 
the same QQq potential in the region of R < 0.8 fm. For 
example, in the M q ~ 1 GeV case, the effective string 
tension is 0.75 ± 0.08 GeV/fm. In order to compare with 
the lattice QCD result, we calculate in the present po- 
tential model with the same condition, M q = 1 GeV and 
R < 0.8 fm, and then find a cS ~ 0.76 GeV/fm. There- 
fore, the calculation of our simple potential model almost 
reproduces the result of lattice QCD. 



B. The reduction of the effective string tension 

To understand the reduction mechanism of a e g be- 
tween the two heavy quarks, we compare the definition 
of the string tension with that of the effective string ten- 
sion. The string tension is the proportionality coefficient 
of the color flux-tube length L m { n in the confinement po- 
tential. Compared with this, the effective string tension 
is defined as the effective confinement force between two 
quarks in baryons in terms of the inter-two-quark dis- 
tance R . The schematic figure is depicted in Fig. [4] In 
the QQ case, L m i n = R and the two definitions lead no 
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R R 

FIG. 4: The schematic figure of the flux-tube length L m i n 
and the inter-quark distance R. The two definitions are the 
same, L m i n = R, in the QQ case (left). These are different, 
imin 7^ R, in the 3Q or QQ<j case (right). 



difference, but in the 3Q or QQq case, these are different 
from each other, L min ^ R. The QQq flux-tube length is 
related to R through complicated light-quark dynamics, 
and then the QQq effective confinement potential is not 
necessarily a linear function of R. Our results show that 
the effective confinement potential is almost proportional 
to R, at least in the region of R < 1.2 fm, with the pro- 
portionality coefficient <7 c ff, which is reduced compared 
to the string tension of mesons and the static cases. We 
conjecture that the reduction of a e s originates from this 
difference between the inter-heavy-quark distance R and 
the flux-tube length L m i n - 

The relation between R and the QQq flux-tube length 
can be easily calculated in our potential model. The ex- 
pectation value of the QQq flux-tube length (L m i n ) is 
shown in Table Q] and Fig. [5l As shown in Fig. [5] (solid 
line), (L m in) is well described as (L m i n ) ~ bo + b\R in 
the range of R < 1.2 fm, and its best-fit parameters 
are bo ~ 0.61 fm and b± ~ 0.81. Apart from the light- 
quark kinetic energy and the Coulomb contribution, wc 
can explain the reduction of a e g by the functional form 
of (Lmin), that is, 

(&3QL m in) ~ (bl<7 3 Q)R + bo<JZQ 

~ fToff-R + const. (16) 

Thus (T e ff is reduced as cr e ff ~ b\a^Q, and &i(< 1) means 
the reduction rate of <7 e ff- 

From the different M q results, we see the ./In- 
dependence of the system. Intuitively, in the M q — > oo 
limit, the QQq system becomes the static 3Q system and 
then CTeff approaches to 03Q. It is expected that, when M q 
increases, the reduction is weakened and <7 e ff increases. 
Similarly, since quarks are not spatially extended at all in 
the 3Q case, the light-quark spatial distribution becomes 
compact for large M q . We can check these behaviors in 
Table U and HB 

Now we investigate the contribution of the Coulomb 
terms in Eq. ([3]). We artificially set A 3 q — 0, and cal- 
culate in the same scheme as above. We fit this "No- 
Coulomb" QQq potential with a functional form 

V$S q (R) = <r? s c R + C? s c . (17) 




0.5 - 



1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 

R [fm] 

FIG. 5: The relation between (L min ) and R with M q = 330 
MeV. The symbol + denotes the result from the variational 
calculation. The solid line is the best-fit function of 60 + b\R 
with the fit range of R < 1.2 fm, and the dashed line is that 
of bo + biR + b 2 R 2 with R < 2.4 fm. 



TABLE IV: The "No-Coulomb" (A 3Q = 0) QQq potential 
results with the 128 x 128-mesh. cr^j and C^f are the best- 
fit parameters in Eq. (|17[1 with the fit range R < 1.2 fm. 



M q [GeV] 


Az [fm] 


a c N ff c [GeV/fm] 


C P N ff c [GeV] 


0.33 


0.10 


0.69 


1.19 




0.12 


0.69 


1.19 


0.50 


0.08 


0.70 


1.25 




0.10 


0.70 


1.25 



The best-fit parameters are shown in Table HVl Also in 
this case, the effective string tension afg is reduced, and 
roughly equals to a e g. Therefore, the essential reason for 
the reduction of cr e ff is the geometrical difference between 
Lmin and R rather than the Coulomb contribution. 

In Fig. [51 we separate the QQq potential with M q = 
330 MeV into the contributions from the expectation 
value of the light-quark kinetic term, the 3Q confinement 
potential term, and the Coulomb term. The resulting 
QQq potential is given as V QQq (R) = M q + (-5^7^ + 

03QL m in — J2i<j We can confirm that the confine- 

ment part of the QQq potential originates mainly from 
the 3Q confinement potential contribution. 

C. The detail of the effective string tension 

We note that the QQq potential is not perfectly fitted 
with Eq. (fl4|) . For example, in the large R limit, R will 
approach to the flux-tube length (L m i n ), so that a e g in- 
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FIG. 6: The component of the QQq potential with M q = 
330 MeV. The solid curve is the best-fit function VQQ q (R) 
of Eq. g4j), and it is given as V QQq (R) = M q + (f /2M q + 



creases and approaches to o~ 3 q ■ Thus a c ff should depend 
on R. To investigate the more accurate i?-dependence, 
we enlarge the region of R to 2.4 fm. As mentioned above, 
(£min) is determined by a nontrivial light-quark dynam- 
ics. Then, in general, it is a function with higher-power 
terms of R. 



(imin) 

If we fit (L r 



= 6 + hR + b 2 R 2 + b 3 R 3 



with 



(18) 



M q = 330 MeV and R < 2.4 fm 



up to the quadratic term with 
we find that the best-fit 
parameters are b ~ 0.62 fm, b\ ~ 0.77, and b 2 — 0.037 
fm -1 . Its best-fit function is shown in Fig. [5] (dashed 
line). In this region, the linear function deviates from 
the result, and this quadratic function seems suitable for 
fitting. In this case, if the light-quark kinetic energy and 
the Coulomb contribution is ignored, the effective string 
tension is roughly written as 



Ceff(-R) 



_ dV QQq {R) 



dR 



bia 3 g + 2b 2 a 3Q R, (19) 



where "IR" means the infrared region. Thus the effec- 
tive string tension actually depends on R. To investigate 
the independent effective string tension with the whole 
interaction of Eq. ([3]) , we fit the QQq potential with 



A 



off 



V QQq (R) = (&! + b 2 R)a 3Q R - -|i 



a 



cir 



(20) 



in the range of R < 2.4 fm. The results are summarized 
in Table fVl From the fact that b 2 is positive, a e s(R) is 
an increasing function of R, as is expected. Of course, 
this argument is restricted in the present range of R. For 
larger i?, higher-power terms are to be included. 



TABLE V: The best-fit parameters of Eq. (J20J with the fit 
range of R < 2.4 fm. The calculation is done with M q = 330 
MeV and the 128 x 128- mesh. 



Az [fm] 


61 


b 2 [fm" 1 ] 


A eS 


C eff [GeV] 


0.10 


0.78 


0.031 


0.12 


1.08 


0.12 


0.78 


0.033 


0.12 


1.08 



D. Comments on the relativistic corrections 

The non-relativistic quark model has achieved consid- 
erable success in explaining low-energy hadron proper- 
ties. However, the relativistic treatment for the light 
quark can give rise to the quantitative correction to some 
properties. We give a brief comment on the relativistic 
correction to our result. 

In the relativistic formalism for the light quark, the 
Hamiltonian |(3J) is modified as 



H 



^/M2+p 2 + y(f 1 ,r 2 ,r 3 ), (21) 

where p is the light-quark momentum. The interaction 
term V(fi, r 2 ,f 3 ) is again the general three-quark inter- 
action in QQq systems. 

One of the relativistic corrections is the higher-order 
kinetic term from (AI 2 q + p 2 ) 1 / 2 - (M q +p 2 /2M q ). If the 
higher-order terms do not change the i?-dependence of 
the QQq potential, the non-relativistic formalism works 
well. As shown in Fig. the leading-order kinetic term 
p 2 /2M q is almost i?-independent, compared to the typ- 
ical i?-dependence of the energy, namely <t 3 q — 0.89 
GeV/fm. We expect that the higher-order kinetic terms 
also give rise to small i?-dependent contributions, at least 
in the large M q case. 

Another possible relativistic correction is modification 
to V(fi,f 2 ,r 3 ). While we adopt the static 3Q potential 
as the three-quark interaction in QQq systems, it can 
be changed by the light-quark relativistic effect. Several 
works have dealt the inter-quark potential in mesons and 
baryons with finite-mass quarks [l5|, [TtJ ■ The finite-mass- 
quark effect on the potential is rather complicated, and 
nontrivial Mq-dependence can occur. 

For the precise statement and quantitative accuracy, 
the relativistic treatment is desired, particularly for the 
realistic light-quark mass case, i.e., M q ~ 330 MeV. 



IV. SUMMARY AND CONCLUSION 

From the above considerations, we can conclude that, 
when finite-mass quarks exist in baryons, the effective 
string tension between the other two quarks is reduced 
compared to the static case. The essential reason for 
the reduction is the geometrical difference between the 
inter-quark distance R and the flux-tube length (L m i n ). 
As can be seen from Fig. 21 this geometrical difference al- 
ways exists for more than two-quark system, regardless of 
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whether the heavy-quark mass is infinite or finite. Then 
this is expected to hold not only for QQq systems but 
also for ordinary baryons, such as nucleons. Of course, 
the quantitative difference would exist between QQq sys- 
tems and ordinary baryons. In addition, also in multi- 
quark systems, the effective string tension can be changed 
due to the existence of other light quarks. In such sys- 
tems, the effective string tension between two finite-mass 
quarks would include more complicated finite-mass ef- 
fects. 

In summary, we have studied the QQq potential with 
the non-relativistic potential model, and have investi- 
gated the finite-mass-quark motional effect on the inter- 
heavy-quark potential in baryons. We have found the 
significant reduction of the effective string tension a e s be- 



tween the two heavy quarks in QQq systems, compared 
to string tension u^q in the static 3Q case. The effec- 
tive string tension a e s depends on the light-quark mass 
M q , and slightly depends on R. The finite-mass-quark 
existence reduces the effective confinement force, which 
is conjectured to be a general property for baryons. 
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